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INTRODUCTION 

ven  with  all  the  high  technology 
weapons  onboard  U.S.  Navy  ships  today, 
the  difference  between  success  and  failure 
often  comes  down  to  our  understanding 
and  knowledge  of  the  environment  in 
which  we  are  operating.  Accurately  pre¬ 
dicting  the  ocean  environment  is  a  criti¬ 
cal  factor  in  using  our  detection  systems 
to  find  a  target  and  in  setting  our  weap¬ 
ons  to  prosecute  a  target  (Gottshall,  1997; 
Chu  et  al.,  1998).  From  the  ocean  tem¬ 
perature  and  salinity,  the  sound  velocity 
profiles  (SVP)  can  be  calculated.  SVPs  are 
a  key  input  used  by  U.S.  Navy  weapons 
programs  to  predict  weapon  performance 
in  the  medium.  The  trick  lies  in  finding 
the  degree  to  which  the  effectiveness  of 
the  weapon  systems  is  tied  to  the  accu¬ 
racy  of  the  ocean  predictions. 

The  U.S.  Navy’s  Meteorological  and 
Oceanographic  (METOC)  community 
currently  uses  three  different  methods  to 
obtain  representative  SVPs  of  the  ocean: 
climatology,  in  situ  measurements,  and 
data  (including  satellite  data)  assimilation. 
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Impact  of  satellite  data  assimilation  on  naval  undersea  capability  is  investigated 
using  ocean  hydrographic  products  without  and  with  satellite  data  assimilation.  The 
former  is  the  Navy’s  Global  Digital  Environmental  Model  (GDEM),  providing  a  monthly 
mean;  the  latter  is  the  Modular  Ocean  Data  Assimilation  System  (MODAS)  providing 
synoptic  analyses  based  upon  satellite  data.  The  two  environmental  datasets  are  taken 
as  the  input  into  the  Weapon  Acoustic  Preset  Program  to  determine  the  suggested 
presets  for  an  Mk  48  torpedo.  The  acoustic  coverage  area  generated  by  the  program 
will  be  used  as  the  metric  to  compare  the  two  sets  of  outputs.  The  output  presets 
were  created  for  two  different  scenarios,  an  anti-surface  warfare  (ASUW)  and  an 
anti-submarine  warfare  (ASW);  and  three  different  depth  bands,  shallow,  mid,  and 
deep.  After  analyzing  the  output,  it  became  clear  that  there  was  a  great  difference  in 
the  presets  for  the  shallow  depth  band,  and  that  as  depth  increased,  the  difference 
between  the  presets  decreased.  Therefore,  the  MODAS  product,  and  in  turn  the  sat¬ 
ellite  data  assimilation,  had  greatest  impact  in  the  shallow  depth  band.  The  ASW 
presets  also  seemed  to  be  slightly  less  sensitive  to  differences  than  did  presets  in 
the  ASUW  scenario. 


The  climatological  data  provides  the  back¬ 
ground  SVP  information  that  might  not 
be  current.  The  Generalized  Digital  En¬ 
vironmental  Model  (GDEM)  is  an  ex¬ 
ample  of  a  climatological  system  that  pro¬ 
vides  long-term  mean  temperature,  salin¬ 
ity,  and  sound  speed  profiles.  The  in  situ 
measurements  from  conductivity-tem¬ 
perature-depth  (CTD)  and  expendable 
bathythermograph  (XBT)  casts  may  give 
accurate  and  timely  information,  but  these 
are  not  likely  to  have  large  spatial  and  tem¬ 
poral  coverage  over  all  regions  where  U.S. 
ships  are  going  to  be  operating.  In  a  data 
assimilation  system,  an  initial  climatology 
or  forecast  is  improved  by  using  satellite 
and  in  situ  data  to  better  estimate  synop¬ 
tic  SVPs.  The  Modular  Ocean  Data  As¬ 
similation  System  (MODAS)  utilizes  sea 
surface  height  (SSH)  and  sea  surface  tem¬ 
perature  (SST)  in  this  way  to  make 
nowcasts  of  the  ocean  environment  (Fox 
et  al„  2002). 

The  value  added  by  satellite  data  as¬ 
similation  for  use  of  undersea  weapon  sys¬ 
tems  can  be  evaluated  using  the  SVP  in¬ 


put  data  from  MODAS  (with  satellite  data 
assimilation)  and  GDEM  (climatology 
without  satellite  data  assimilation).  The 
question  also  arises  of  how  many  altim¬ 
eters  are  necessary  to  generate  an  optimal 
MODAS  field.  Too  few  inputs  could  re¬ 
sult  in  an  inaccurate  MODAS  field,  which 
in  turn  could  lead  to  decreased  weapon 
effectiveness.  There  must  also  be  some 
point  at  which  the  addition  of  another 
altimeter  is  going  to  add  a  negligible  in¬ 
crease  in  effectiveness.  This  superfluous 
altimeter  would  be  simply  a  waste  of 
money  that  could  be  spent  on  a  more  use¬ 
ful  system. 

The  purpose  of  this  study  is  to  quan¬ 
tify  the  advantage  gained  from  the  use  of 
data  from  MODAS  assimilation  of  satel¬ 
lite  observations  rather  than  climatology. 
The  study  will  specifically  cover  the  ben¬ 
efits  of  MODAS  data  over  climatology 
when  using  their  respective  SVPs  to  deter¬ 
mine  torpedo  settings.  These  settings  re¬ 
sult  in  acoustic  coverage  percentages  that 
will  be  used  as  the  metric  to  compare  the 
two  types  of  data. 
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2.  Navy’s  METOC  Models 
and  data 

2.1.  GDEM 

GDEM  is  a  four  dimensional  (latitude, 
longitude,  depth  and  time)  digital  model 
maintained  by  the  Naval  Oceanographic 
Office.  GDEM  was  generated  using  over 
seven  million  temperature  and  salinity  ob¬ 
servations,  most  of  them  drawn  from  the 
Master  Oceanographic  Observation  Data 
Set  (MOODS).  Globally  GDEM  has  a  reso¬ 
lution  of  1/2°  degree.  However,  in  a  few  se¬ 
lect  areas,  higher  resolutions  are  available. 
In  order  to  represent  the  mean  vertical  dis¬ 
tribution  of  temperature  and  salinity  for  grid 
squares,  GDEM  determines  analytical 
curves  to  fit  to  the  individual  profiles  (Teague 
et.  al„  1990;  Chu  et  al„  1997;  1999) 

Before  curves  can  be  fitted  to  the  data, 
quality  control  must  be  implemented  that 
removes  anomalous  features  or  bad  obser¬ 
vations.  The  data  is  checked  for  proper  range 
and  static  stability,  and  it  is  checked  to  en¬ 
sure  that  it  has  not  been  misplaced  in  loca¬ 
tion  or  season.  Once  the  data  has  been  in¬ 
spected  for  quality,  curves  are  fitted  to  the 
data.  From  the  mathematical  expressions 
that  represent  the  curves,  coefficients  are 
determined.  It  is  these  coefficients  that  will 
be  averaged.  It  can  be  shown  that  the  coeffi¬ 
cients  resulting  from  averaged  data  are  not 
the  same  as  the  averaged  coefficients  of  the 
data.  In  order  to  minimize  the  number  of 
coefficients  necessary  to  generate  smooth 
curves,  different  families  of  curves  are  used 
for  different  depth  ranges.  This  necessitates 
the  careful  selection  of  matching  conditions 
in  order  to  ensure  that  no  discontinuities  in 
the  vertical  gradients  occur.  Separate  com¬ 
putation  of  temperature  and  salinity  allow 
the  results  to  be  checked  against  each  other 
to  ensure  stable  densities. 

2.2.  MODAS 

MODAS  is  a  collection  of  over  100 
FORTRAN  programs  and  UNIX  scripts 
that  can  be  combined  to  generate  a  number 
of  different  products  (Fox  et  al.,  2002).  A 
few  examples  of  MODAS  programs  include 
data  sorting,  data  cross-validation,  data  as¬ 
similation,  and  profile  extension.  This 


modularity  allows  MODAS  to  be  quickly 
and  easily  modified  to  handle  problems  or 
new  requirements  as  they  arise.  MODAS 
has  varying  degrees  of  resolution  starting  at 
1/2°  in  the  open  ocean  increasing  to  1/4°  in 
coastal  seas  and  increasing  again  to  1/8°  near 
the  coast  (Fox  et  al.,  2002).  To  generate 
nowcasts  and  forecasts,  the  MODAS  sys¬ 
tem  uses  a  relocatable  version  of  the 
Princeton  Ocean  Model  (POM).  To  initial¬ 
ize  the  POM  MODAS  temperature  and 
salinity  grids,  geostrophically  estimated  cur¬ 
rents,  or  extracted  currents  from  other 
POM’s  can  be  used. 

One  of  the  most  important  features  of 
MODAS  is  its  use  of  dynamic  climatology 
(Fox  et  al.,  2002).  Dynamic  climatology  is 
the  incorporation  of  additional  information 
into  the  historical  climatology  in  order  to 
portray  transient  features  that  are  not  repre¬ 
sented  by  the  climatology.  Two  useful  quan¬ 
tities  that  are  easily  gathered  from  satellites 
are  sea  surface  height  (SSH)  and  sea  surface 
temperature  (SST).  While  SST  from  altim¬ 
eters  can  be  used  directly,  the  SSH,  which  is 
measured  as  the  total  height  relative  to  the 
proscribed  mean,  must  be  converted  into  a 
steric  height  anomaly  in  order  to  be  used. 
2D  SST  and  SSH  fields  are  generated  from 
point  observations  through  the  use  of  opti¬ 
mal  interpolation. 

Optimal  interpolation  is  a  process  by  which 
the  interpolated  temperature  or  salinity 
anomaly  is  determined  as  the  linear  combi¬ 
nation  of  the  observed  anomalies.  Each  of 
the  anomalies  is  given  a  weight  that  accounts 
for  variation  in  temporal  and  spatial  sam¬ 
pling.  Weights  are  computed  by  minimiz¬ 
ing  the  least  square  difference  between  the 
interpolated  value  and  the  true  value  at  the 
grid  point  and  by  solving  the  equations 
(Gandin,  1965), 

N 

'ZajMij+^2ai=MGi,  (1) 

7=1 

where  a(  are  the  weights,  A  is  the  signal  to 
noise  ratio,  is  the  autocorrelation  be¬ 
tween  locations  i  and  j,  and  fiGj  is  the 
autocorrelation  between  the  grid  point  and  i. 
For  each  grid  node  location  matrix  inversion 
is  used  to  solve  the  system  of  N  equations 


for  the  N  unknown  weights.  The  other  pa¬ 
rameters  are  computed  using  the  first  guess 
field,  MOODS  profiles,  and  climatology. 
Using  this  process  any  new  observation  can 
be  interpolated  into  the  appropriate 
MODAS  grid  node. 

The  first  guess  field,  the  prior  day’s  2D 
SST  field,  or  the  weighted  average  of  35  days 
of  altimeter  data  respectively,  is  subtracted 
from  the  new  observations,  and  the  result¬ 
ing  deviations  are  interpolated  to  produce  a 
field  of  deviation.  This  is  added  to  the  first 
guess  field  to  generate  the  new  2D  field.  For 
the  first  iteration  of  the  optimal  interpola¬ 
tion,  climatology  is  used  for  SST  and  the 
SSH  measurement  is  assumed  to  have  a  zero 
deviation.  This  means  that  until  the  field 
deviates  from  the  climatology,  the  extra  data 
has  added  no  value  and  MODAS  reverts  to 
climatology. 

Once  the  data  is  in  a  useful  form, 
MODAS  begins  with  the  climatology  pro¬ 
file  and  then  correlates  variations  in  the  SSH 
and  SST  to  variations  in  the  subsurface  tem¬ 
perature.  The  regression  relationships  used 
here  were  constructed  by  performing  a  least- 
squares  regression  analysis  on  archived  tem¬ 
perature  and  salinity  profiles.  This  is  a  three 
step  process  starting  with  the  computation 
of  regional  empirical  orthogonal  functions 
from  the  historical  temperature  and  salinity 
profiles.  The  second  step  is  to  express  the 
profiles  in  terms  of  an  empirical  orthogonal 
function  series  expansion.  The  final  step  is 
to  perform  regression  analysis  on  the  profile 
amplitudes  for  each  mode,  truncating  the 
series  after  three  terms.  This  is  possible  be¬ 
cause  of  the  compactness  of  the  empirical 
orthogonal  function  representation. 

Once  the  subsurface  temperatures  have 
been  revised,  MODAS  adjusts  the  subsur¬ 
face  salinity  profile  using  the  relationship 
between  temperature  and  salinity.  This 
new  profile  is  referred  to  as  a  synthetic 
profile.  Synthetic  profiles  only  utilize  these 
regression  relationships  down  to  a  depth 
of  1500  m  due  to  the  decreasing  reliabil¬ 
ity  of  the  relationships  at  depth  (Fox  et 
al.,  2002). 

MODAS  is  also  able  to  include  mea¬ 
surements  from  in  situ  CTDs  and  XBTs.  The 
first  guess  field  is  the  field  generated  by  the 
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dynamic  climatology,  and  the  in  situ  pro¬ 
files  are  subtracted  from  it  to  get  residuals. 
Optimal  interpolation  is  once  again  used  to 
update  the  temperature  field  and  from  the 
temperature  field  the  salinity  field  can  be 
generated.  This  salinity  field  then  serves  as  a 
first  guess  field  for  the  inclusion  of  the  salin¬ 
ity  profiles  (Fox  et  al.,  2002). 

2.3  Satellite  Altimetry  Data 
Assimilated  into  MODAS 

The  Navy  currently  uses  satellite  al¬ 
timeters  and  inferred  data  to  assess  the 
ocean  environment  for  the  naval  opera¬ 
tions.  Of  primary  interest  is  mesoscale 
variability.  Meandering  fronts  and  eddies 
can  significantly  change  the  temperature 
and  salinity  structure  of  the  ocean.  This 
importance  is  clearly  seen  in  sonar  depen¬ 
dent  operations  such  as  anti-submarine 
warfare  (ASW).  Sonar  range  can  be  greatly 
helped  or  hindered  by  the  acoustic  envi¬ 
ronment  created  by  the  salinity,  tempera¬ 
ture,  and  density.  Altimeters  also  provide 
the  SSH  measurements  that  MODAS  uses 
in  its  optimal  interpolation. 

While  monitoring  mesoscale  variability 
is  of  prime  importance  to  the  Navy,  an 
emerging  secondary  role  for  Navy  altimeters 
is  monitoring  continental  shelf  and  coastal 
zones.  As  the  Navy  conducts  more  and  more 
operations  in  littoral  waters,  the  ability  to 
predict  near-shore  parameters  will  have  in¬ 
creasing  importance.  Altimeter  data  can  be 
used  to  get  up-to-date  information  on  rap¬ 
idly  changing  near-shore  characteristics  such 
as  tides  and  wave  height  (Jacobs  et  al.,  2002). 
These  are  important  issues  for  anyone  deal¬ 
ing  with  mine  detection,  beach  operations, 
or  ship  routing. 

Altimeters  have  also  been  used  to  mea¬ 
sure  the  flow  through  important  straits,  such 
as  the  Tsushima  Strait,  and  to  measure  large- 
scale  circulation.  The  first  of  these  helps  re¬ 
searchers  and  modelers  to  develop  con¬ 
straints  on  local  numerical  models.  Large- 
scale  circulation  measurements  can  also  help 
in  the  development  of  models  by  aiding  in 
error  correction.  They  also  help  explain  the 
local  environment  that  is  often  affected  by 
not  just  local  forcing,  but  large-scale  circu¬ 
lation  variations  as  well. 


Satellite  altimeters  can  provide  a  great 
variety  of  data,  but  no  single  altimeter  can 
provide  measurements  on  all  desired  time 
and  length  scales.  Different  parameters  must 
be  sampled  at  different  frequencies  if  they 
are  going  to  be  of  any  use.  For  instance,  sea 
surface  height  must  be  sampled  every  48 
hours  while  wave  height  must  be  sampled 
every  three  hours.  While  different  ocean  fea¬ 
tures  all  have  different  time  and  spatial  scales, 
only  the  requirements  for  observation  of 
mesoscale  features  are  presented  here  as  an 
example  (Jacobs  et.  al.,  1999). 

In  order  for  an  altimeter  to  efficiently 
and  accurately  sample  mesoscale  features, 
there  are  several  requirements  placed  on  its 
accuracy,  orbit,  and  repeat  period.  A  satel¬ 
lite  altimeter  must  produce  measurements 
that  are  accurate  to  within  5  cm,  or  the  er¬ 
rors  that  propagate  down  into  the  tempera¬ 
ture  and  salinity  calculations  will  be  unac¬ 
ceptable.  With  an  error  of  only  5  cm,  the 
error  in  the  temperature  calculation  can  be 
1-2°  C.  Satellites  should  also  have  an  ex¬ 
act  repeat  orbit  to  maximize  the  usefulness 
of  the  data  collected.  Without  an  exact  re¬ 
peat  orbit,  the  only  way  to  get  differences 
in  sea  surface  heights  is  to  use  only  the  data 
from  points  where  the  satellite  crosses  the 
track  of  another  altimeter  or  itself.  An  ex¬ 
act  orbit  is  considered  to  be  a  1  km  wide 
swath  of  a  predefined  ground  track.  Fi¬ 
nally,  the  period  of  a  single  satellite  should 
be  greater  than  the  typical  20  day  time  scale 
of  a  mesoscale  feature.  If  two  satellites  are 
used,  then  they  should  be  spaced  so  that  a 
point  on  the  ground  is  not  sampled  more 
than  once  in  a  20  day  period  (Jacobs  et. 
al.,  1999). 

As  described  earlier,  systems  such  as 
MODAS  rely  heavily  on  the  information 
provided  by  these  satellites.  MODAS  uses 
interpolation  to  estimate  SSFf  at  points 
that  the  satellite  did  not  cover.  If  the 
ground  track  spacing  is  too  coarse  then 
the  optimal  interpolation  scheme  of 
MODAS  will  begin  introducing  errors 
into  the  fields  between  the  tracks.  It  is 
important  that  the  satellites  be  properly 
set  up  so  that  a  maximum  amount  of  in¬ 
formation  can  be  gathered  with  a  mini¬ 
mum  amount  of  error. 


3. Navy’s  Weapon  Acoustic  Preset 

A  Weapon  Acoustic  Preset  Program 
(WAPP)  is  used  to  get  automated,  interac¬ 
tive  means  of  generating  Mk  48  and  Mk  48 
Advanced  CAPability  (ADCAP)  acoustic 
presets  and  visualizing  torpedo  performance. 
It  combines  the  Mk  48  Acoustic  Preset  Pro¬ 
gram  (M48APP)  and  the  Mk  48  ADCAP 
Acoustic  Preset  Program  (MAAPP)  into  a 
single  integrated  package.  The  Royal  Aus¬ 
tralian  Navy  as  a  part  of  the  Collins  Class 
Augmentation  System  (CCAS)  also  uses  the 
M48APP,  and  the  Royal  Canadian  Navy  has 
changed  the  M48APP  for  Java.  The  program 
is  based  around  a  graphical  user  interface 
that  allows  the  user  to  enter  the  environmen¬ 
tal,  tactical,  target,  and  weapon  data.  With 
these  user  specified  parameters,  the  program 
then  performs  a  series  of  computations  to 
generate  accurate  acoustic  performance  pre¬ 
dictions.  The  output  includes  a  ranked  list- 
set  of  search  depth,  pitch  angle,  LD,  and 
effectiveness  values,  an  acoustic  ray  trace,  and 
a  signal  excess  map  (Cwalina,  2002,  per¬ 
sonal  communication). 

The  Environmental  Data  Entry  Mod¬ 
ule  (EDE)  is  a  simple  Graphic  User  Inter¬ 
face  (GUI)  that  allows  the  user  to  enter  a 
variety  of  environmental  parameters  (Fig.  1). 
The  sea  surface  fields  allow  the  user  to  specify 
wind  speed,  wave  height,  and  sea  state  based 
on  either  the  World  Meteorological  or  Beau¬ 
fort  scale  conventions.  The  three  fields  are 
coupled  so  that  an  entry  into  one  field  will 
bring  up  the  appropriate  default  values  for 
the  others.  The  bottom  condition  field  al¬ 
lows  the  user  to  specify  the  bottom  depth 
and  to  choose  the  bottom  type  from  a  list  of 
possibilities.  The  bottom  of  the  GUI  is  de¬ 
voted  to  the  water  column  characteristics  and 
a  sound  speed  profile.  The  temperature, 
sound  speed,  and  depth  are  all  in  the  appro¬ 
priate  English  units.  The  volume  scattering 
strength  (VSS)  is  in  dB.  The  additional  fields 
include  the  latitude,  longitude,  the  profile 
name,  and  the  table  group  identifiers. 

Once  the  environmental  parameters 
have  been  entered,  the  user  can  move  on  to 
the  Acoustic  Module  Preset  Display.  This 
GUI  allows  the  user  to  specify  a  number  of 
parameters  about  the  weapon,  the  target,  and 
the  way  the  weapon  should  search  (Fig.  2). 
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FIGURE1 

Environmental  data  entry. 


FIGURE  2 

Acoustic  preset  module  display. 


The  list-set  on  the  right  side  of  the  GUI  dis¬ 
plays  a  series  of  search  depths,  pitch  angles, 
laminar  distances,  and  effectiveness  values. 
The  effectiveness  values  for  the  various  pre¬ 
sets  are  based  on  expected  signal  excess  and 
ray  trace  computations.  Both  plots  can  be 
viewed  from  a  pull-down  menu.  These  pro¬ 
vide  a  visual  representation  of  the  acoustic 
performance  of  the  Mk  48. 

In  addition  to  automatically  computing 
the  most  effective  preset  combination  for  a 
given  set  of  environmental  parameters,  the 
program  also  allows  the  user  to  manually 
examine  the  effectiveness  of  any  allowable 
preset  combination  via  the  signal  excess  and 
ray  trace  plots.  The  program  also  allows  the 
user  to  save  the  tactical  preset  list  and  the 
accompanying  environmental  data.  The  data 
are  stored  locally  in  the  weapon  module  and 
can  be  recalled  later  or  transferred  via  a  net¬ 
work  to  the  combat  control  system. 


4.  Statistical  Analysis 

4.1.  Input  and  Output  Difference 

The  difference  between  the  two  sets  of 
input  (GDEM  and  MODAS)  I jf  input ,  or 
between  the  two  sets  of  output  weapon  pre¬ 
set  data  (running  using  GDEM  and 
MODAS)  1 \foutput, 

Aiy(r ,t)  =  -  y/G(r ,t) ,  (2) 

represents  the  ocean  data  update  using  sat¬ 
ellite  and  in  situ  observations  (input)  and 
the  effect  of  using  satellite  and  in  situ  obser¬ 
vations  on  the  weapon  preset  (output).  Here 
vv  and  y/G  are  the  variables  (either  input 
or  output)  using  GDEM  and  MODAS,  re¬ 
spectively.  We  may  take  the  probability  his¬ 
tograms  of  and  to  show  the  difference  be¬ 
tween  the  statistical  characteristics. 

4.2.  Root  IVlean  Square  Difference 

GDEM  and  MODAS  have  different 
grid  spacing:  1/2°  1/2°  in  GDEM  and 
1  / 1 2°  1  / 1 2°  in  MODAS.  For  a  GDEM  cell, 
one  data  is  available  for  GDEM  and  36  data 
for  MODAS.  The  root-mean-square  differ¬ 
ence  (RMSD), 


14  Marine  Technology  Society  Journal 


RMSD  : 
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N 


,  (3) 


is  commonly  used  to  represent  the  differ¬ 
ence  in  the  input  and  output  data.  Here,  N 
(=36)  is  the  total  MODAS  data  number  in 
a  GDEM  cell.  The  RMSD  can  be  computed 
for  either  the  input  data  to  the  weapon  pre¬ 
set  model  such  as  the  temperature,  salinity, 
or  sound  speed,  or  it  can  be  computed  for 
the  output  data  such  as  nondimensional  de¬ 
tection  area. 


Due  to  the  differing  resolutions  of 
GDEM  and  MODAS,  this  area  provided 
117  GDEM  profiles  and  1633  MODAS 
profiles  (Fig.  4).  Each  profile  was  simply  a 
text  file  that  consisted  of  a  header  row  and 
columns  of  data.  The  header  row  contained 
the  number  of  depths  the  profile  covered, 
the  file’s  name  and  the  latitude  and  longi¬ 
tude  of  the  profile. 

FIGURE  4 

GDEM  and  MODAS  data  points 


The  columns  corresponded  to  depth  in 
feet,  the  temperature  in  degrees  Fahrenheit, 
the  sound  speed  velocity  in  feet  per  second, 
a  volume  backscatter  value,  and  salinity  in 
PSU.  Despite  the  common  use  of  Interna¬ 
tional  units  in  scientific  experiments,  it  was 
necessary  for  the  profiles  to  be  set  up  in  the 
appropriate  English  units.  The  Weapon 
Acoustic  Preset  Program  (WAPP),  the  pro¬ 
gram  used  to  generate  the  presets  from  the 
profiles,  requires  inputs  to  be  in  English  units. 


5.  Comparison  between 
gdem  and  modas  in  the 
gulf  stream  region 

5.1.  Data 

In  order  to  make  a  meaningful  compari¬ 
son  of  MODAS  and  GDEM  data,  a  suffi- 
ciendy  large  data  set  had  to  be  obtained.  The 
Area  of  Interest  (AOI)  also  needed  to  be  an 
area  where  the  ocean  environment  fluctuated 
on  a  fairly  short  time  scale.  The  GDEM  data 
in  March  and  MODAS  data  on  March  15, 
2001  was  obtained  for  the  area  off  the  North 
American  coast  corresponding  to  40°-35°  N 
latitude  and  75°-70°  W  longitude  (Fig.  3). 

FIGURE  3 

Area  of  interest. 


5.2.  Difference  between  GDEM 
and  MODAS 

While  GDEM  gives  the  climatological 
background  ocean  environment  at  a  given 
place,  MODAS  is  expected  to  provide  more 
current  and  synoptic  interpretations  of  the 
environment.  The  amount  of  accuracy 
MODAS  adds  is  in  proportion  to  the  scale 
on  which  ocean  parameters  vary.  For  areas 
such  as  the  Gulf  Stream,  where  environmen¬ 
tal  factors  are  known  to  vary  rapidly  on  a 
relatively  small  time  scale,  it  is  expected  that 
there  would  be  at  least  a  few  areas  where  the 
two  data  sets  differ.  It  is  these  areas  that  are 
of  particular  interest,  since  the  difference  in 
the  weapon  presets  should  be  greatest. 

On  the  surface,  the  GDEM  data  pro¬ 
vided  a  view  of  the  temperature  distribution 
that  consisted  of  smooth,  uniformly  spaced 
lines  of  constant  temperature  that  were  con¬ 
sistent  with  the  overall  flow  of  the  region 
(Fig.  5).  The  cool  water  on  the  shelf  gradu¬ 
ally  gives  way  to  the  warm  water  flowing 
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FIGURE  5 

March  surface  temperature  and  salinity  distribution  from  GDEM. 
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north  along  the  Gulf  Stream.  The  GDEM 
generated  surface  salinity  distribution  is  ex¬ 
tremely  similar  to  the  surface  temperature 
distribution  and  is  consistent  with  the  Gulf 
Stream  region.  Fresher  water  lies  inland  and 
the  salinity  increases  with  distance  from  the 
shore.  The  only  variation  is  in  the  north¬ 
eastern  section,  where  there  is  a  slight  intru¬ 
sion  of  the  salty  offshore  water. 

As  expected,  the  GDEM  and  MODAS 
distributions  are,  overall,  fairly  similar  in  both 
their  range  of  values  and  overall  distribution. 
They  are  similar  to  each  other  in  shape,  and 
both  show  areas  of  cool  fresh  water  near  the 
coast  and  areas  of  warm  salty  water  lying  off¬ 
shore.  There  are,  however,  a  few  differences, 
with  the  intrusion  of  warm  salty  water  in 
the  northeastern  section  of  the  MODAS  fig¬ 
ure  being  the  most  notable.  There  is  also  an 
area  of  high  temperature  in  the  lower  right 
corner  of  the  MODAS  figure  that  does  not 
show  up  in  the  GDEM  figure.  In  general, 
the  MODAS  figure  shows  a  sharper  front 
with  the  water  increasing  in  temperature  and 
salinity  much  more  rapidly  as  the  distance 
from  the  coast  increases  (Fig.  6) .  The  GDEM 
figure  shows  a  gradual  increase  in  tempera¬ 
ture  and  salinity  starting  in  the  top  left  cor¬ 
ner  and  continuing  almost  entirely  down  to 
the  lower  right  corner.  The  MODAS  figure 
shows  the  water  reaching  maximum  tem¬ 
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perature  and  salinity  quickly  and  then  stay¬ 
ing  constant  to  the  lower  right  corner. 

INSERT  FIGURE  6  NEAR  HERE 

While  the  GDEM  and  MODAS  data 
offer  similar  ranges  of  temperatures,  salini¬ 
ties,  and  sound  speeds  at  the  surface,  the 
distribution  of  the  values  is  quite  different. 
The  histograms  in  Figure  7  reveal  that  while 
the  temperature  values  reported  by  both  data 
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sets  are  similar,  the  MODAS  data  has  a 
higher  proportion  of  profiles  located  in  the 
6°-7°  C  range.  The  difference  in  the  salinity 
graphs  is  even  more  drastic  with  the  bulk  of 
the  GDEM  values  located  in  the  middle  of 
the  range  and  the  MODAS  values  split  be¬ 
tween  the  high  and  low  ends  of  the  range. 
The  sound  speed  graph  indicates  that 
MODAS  typically  reports  higher  sound 
speeds  than  does  the  GDEM  data.  This  is 
not  too  surprising,  since  sound  speed  in  the 
upper  water  column  tends  to  be  tied  closely 
to  temperature,  and  the  MODAS  data  in¬ 
dicates  warmer  water  than  the  GDEM  data. 

Increasing  depth  to  50  m  and  then  100 
m,  it  is  clear  to  see  that,  for  temperature,  the 
distribution  of  the  values  over  the  range  for 
both  sets  of  data  is  quite  similar.  There  is 
still  a  slight  preference  in  the  MODAS 
graphs  towards  higher  temperatures,  but  it 
is  not  as  drastic  as  is  seen  on  the  surface. 
Salinity  is  much  the  same,  with  the  differ¬ 
ence  in  shapes  of  the  two  figures  more  a  fac¬ 
tor  of  the  small  number  of  GDEM  profiles 
as  compared  with  the  number  of  MODAS 
profiles.  Sound  speed  is  the  only  area  where 
the  two  data  sets  continue  to  diverge.  From 
the  50  m  and  100  m  sound  speed  figures  it 
is  clear  that,  with  depth,  the  MODAS  data 
indicates  increasing  sound  speed  and  the 
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FIGURE  6 

MODAS  generated  surface  temperature  and  salinity  distribution  on  March  15,  2001. 
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GDEM  data  predicts  some  sort  of  sound 
speed  minimum  at  depth.  This  is  causing 
the  peak  on  the  MODAS  graph  and  the  peak 
on  the  GDEM  graph  to  move  away  from 
each  other  as  depth  increases. 

By  2000  m  the  temperature  and  salinity 
histograms  for  the  two  data  sets  are  virtually 
identical.  At  this  point  any  perceived  differ¬ 
ence  in  the  two  is  solely  a  factor  of  the  dif¬ 
ference  in  the  number  of  profiles  between 
the  two  data  sets.  For  the  sound  speed  fig¬ 
ures  this  is  the  point  of  maximum  separa¬ 
tion.  The  GDEM  data  indicates  low  sound 
speeds  representative  of  a  deep  sound  chan¬ 
nel,  whereas  the  MODAS  data  indicates  that 
the  sound  speed  has  increased  to  this  point. 
After  this  point  the  GDEM  values  begin  ris¬ 
ing  again  to  match  the  MODAS  data. 

While  the  distribution  of  the  values  over 
the  range  is  a  useful  tool  in  examining  the 
inputs,  it  is  the  difference  between  the  in¬ 
puts  that  is  of  real  importance.  Figure  8 
shows  the  RMS  difference  of  the  inputs. 
From  the  surface  temperature  figure  in  Fig¬ 
ure  8,  the  RMS  difference  of  temperature 
peaks  out  in  the  lower  left  corner  of  the  AOI 
at  about  2°  C.  Besides  the  peak,  the  other 
significant  area  is  the  ridge  starting  in  the 
lower  left  corner  and  running  to  the  middle 
top  of  the  figure.  This  corresponds  to  a  nar¬ 
row  region  where  the  GDEM  distribution 
warmed  slower  than  the  MODAS  distribu¬ 
tion  moving  from  the  coast  out  to  sea.  The 
warm  water  intrusion  is  represented  by  the 
gradual  increase  in  height  of  the  ridge.  The 
salinity  difference  at  the  surface  is  nearly  zero 
for  most  of  the  AOI  and  reaches  its  maxi¬ 
mum  value  of  4.5  PSU  along  the  top  of  the 
region.  The  derived  sound  speed  RMS  dif¬ 
ference,  as  expected,  is  smallest  far  from  the 
coast  where  the  difference  in  temperature 
and  salinity  is  smallest  and  increases  towards 
the  coast. 

As  depth  increases,  the  RMS  difference 
in  temperature  and  sound  speed  changes 
slowly,  but  the  difference  in  salinity  drops  off 
quickly.  Neither  the  temperature  nor  sound 
speed  difference  change  significantly,  but  by 
100  meters  the  RMS  difference  for  salinity 
has  gone  down  to  values  of  less  than  .8  PSU. 
From  100  meters  down,  the  temperature  dif¬ 
ference  begins  to  decrease  slowly,  and  by  2000 


FIGURE  7 

Comparison  between  GDEM  and  MODAS  temperature  histograms 
at  (a)  the  surface,  (b)  100  m  depth,  and  (c)  2000  m  depth. 
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FIGURE  8 

Horizontal  dependence  of  RMSD  at  the  surface  between  GDEM  and  MODAS  for 
(a)  temperature,  (b)  sound  speed,  and  (c)  salinity. 
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meters  the  RMS  difference  for  both  tempera¬ 
ture  and  salinity  has  dropped  to  negligible 
levels  for  most  of  the  AOI.  This  is  expected 
since  MODAS  reverts  to  climatology  at 
depth.  Except  for  the  profiles  in  the  north¬ 
western  corner  of  the  AOI  that  did  not  run 
as  deep  as  the  other  profiles  farther  from  the 
coast,  all  the  RMS  difference  vs  depth  pro¬ 
files  were  remarkably  similar.  All  of  the  tem¬ 
perature  differences  showed  either  a  gradual 
decrease  in  the  difference  down  to  about  1000 
meters  or  a  slight  increase  in  the  difference 
immediately  followed  by  a  gradual  decrease 
in  the  difference  down  to  1000  meters.  At 
about  1000  meters  the  temperature  differ¬ 
ences  all  rapidly  dropped  to  near  zero. 

The  sound  speed  profiles  all  show  the 
difference  increasing  down  to  a  maximum 
value  of  60  m/s  at  around  2000  meters.  Af¬ 
ter  that  the  RMS  difference  drops  off,  and 
approaches  zero  by  3000  meters.  The  cause 
of  the  maximum  at  2000  meters  is  lack  of  a 
deep  sound  channel  according  to  the 
MODAS  data.  The  MODAS  profiles  al¬ 
most  all  have  the  sound  speed  steadily  in¬ 
creasing  down  to  the  maximum  depth 
whereas  climatology  indicates  a  sound  speed 
minimum  at  2000  m.  While  there  is  some 
variation  in  how  quickly  the  salinity  differ¬ 
ences  drop  to  near  zero,  they  are  less  than  1 
PSU  by  200  meters.  Shown  in  Figure  9  is  a 
representative  RMS  difference  profile. 


6.  Comparison  of  weapon 
acoustic  preset  using  gdem 
and  modas 

The  raw  data  was  processed  by  the  Na¬ 
val  Underwater  Warfare  Center  (NUWC) 
Division  Newport.  They  received  the  input 
profiles,  ran  them  through  the  WAPP,  and 
generated  the  output.  Percentage  coverage 
was  calculated  based  on  both  surface 
(ASUW)  and  submarine  (ASW)  scenarios. 
The  submarine  scenario  is  a  low  Doppler 
scenario  consistent  with  diesel  submarine 
operations.  The  coverage  percentages  repre¬ 
sent  coverage  in  the  target  depth  band,  ei¬ 
ther  shallow,  mid,  or  deep.  The  coverage 
percentages  were  also  normalized  over  acous¬ 
tic  modes  to  produce  an  output  that  was 
dimensionless. 
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FIGURE  9 

Depth  dependence  of  RMSD  between  GDEM  and  MODAS  for  (a)  temperature,  (b)  sound  speed,  and  (c)  salinity. 


mean  coverage  for  the  AOI.  The  ASW  sce¬ 
nario  yields  similar  results  except  for  the  fact 
that  the  two  means  were  not  even  statisti¬ 
cally  different.  While  this  seems  to  indicate 
that  the  two  data  sets  are  returning  similar 
results,  there  are  some  important  differences. 
First  are  the  outliers  on  the  GDEM  graphs. 
Values  in  the  high  thirties  to  low  fifties  are 
extremely  rare,  yet  the  GDEM  data  indicate 
that  in  at  least  one  location  for  the  ASUW 
scenario  and  several  for  the  ASW  scenario, 
the  weapon  will  perform  to  this  level.  The 
ASW  scenario  also  has  a  rather  significant 
number  of  GDEM  profiles  that  generate 
below  average  coverage  percentages.  This 
would  indicate  that  GDEM  predicts  that 
coverage  will  vary  greatly  with  location.  In 
comparison  the  MODAS  values  for  both 
scenarios  tend  to  be  very  consistent.  Cover¬ 
age  percentage  varies  little  with  location  due 
to  the  fact  that  most  of  the  profiles  He  within 
a  very  narrow  range.  Overall  GDEM  pre¬ 
dicts  excellent  coverage  some  of  the  time  and 
poor  coverage  the  rest  of  the  time.  MODAS 
data  on  the  other  hand,  indicates  that  cover¬ 
age  percentage  will  not  be  excellent  anywhere 
but  the  expected  values  will  be  uniform  over 
the  whole  shallow  depth  band  region. 


6.1.  Output  Distributions 

The  output  provided  by  NUWC  from 
the  WAPP  runs  consisted  of  twelve  differ¬ 
ent  percentage  coverage  groups,  three  depth 
bands  times  two  scenarios  times  the  two  dif¬ 
ferent  types  of  input  data.  For  the  non-SVP 
derived  WAPP  inputs,  consistent  values  were 
used  throughout  the  runs  to  ensure  that  any 
difference  in  the  outputs  was  a  result  of  dif¬ 
ferences  in  the  GDEM  and  MODAS  data. 
For  each  of  the  groups,  basic  statistics  such 
as  mean,  maximum,  minimum,  and  stan¬ 
dard  deviation  were  computed  and  then  the 
data  were  constructed  into  histograms  to  give 
a  visual  representation  of  how  the  data  are 
distributed. 

In  the  shallow  depth  band  ASUW  sce¬ 
nario  both  MODAS  and  GDEM  yield 
mean  coverage  percentages  that  are  very  close 
to  each  other.  While  statistically  the  means 
are  different,  in  real  world  applications  a  few 
percentage  points  difference  is  negligible  (Fig. 
10).  From  a  user’s  standpoint,  this  means 
that  both  sets  of  data  predict  about  the  same 


FIGURE  10 

Shallow  depth  band  coverage  percentage  distributions:  upper  panels  for  GDEM  and  lower  panels  for  MODAS, 
left  panels  for  ASUW  scenario  and  right  panels  for  ASW  scenario. 
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FIGURE  11 

Medium  depth  band  coverage  percentage  distributions:  upper  panels  for  GDEM  and  lower  panels  for  MODAS, 
left  panels  for  ASUW  scenario  and  right  panels  for  ASW  scenario. 


scenario  has  the  larger  predicted  values,  but 
the  values  in  the  ASW  scenario  are  still  on 
the  upper  end  of  what  is  normal.  The 
MODAS  data  predicts  performance  that  is, 
while  not  particularly  bad,  still  much  more 
pessimistic  than  the  GDEM  predictions. 

For  both  scenarios  the  means  of  the 
GDEM  and  MODAS  derived  predictions 
are  statistically  different  with  the  MODAS 
data  providing  the  smaller  mean  in  both  sce¬ 
narios.  Although  the  dispersion  of  the 
GDEM  data  is  large  in  both  scenarios,  the 
data  is  so  heavily  weighted  towards  the  up¬ 
per  end  that  low  GDEM  coverage  percent¬ 
ages  are  average  values  for  the  MODAS  data 
coverage  percentages.  The  MODAS  data 
coverage  percentages  are  once  again  tightly 
grouped;  the  uniformity  of  the  predicted  cov¬ 
erage  percentages  observed  in  the  two  other 
depth  bands  extends  from  the  surface  down 
to  the  selected  maximum  operating  depth. 


The  mid  depth  band  yielded  results  that 
were  similar  in  distribution  to  the  shallow 
depth  band  (Fig.  11).  Across  both  scenarios 
the  mean  coverage  of  the  GDEM  data  and 
the  mean  coverage  of  the  MODAS  data  are 
statistically  identical.  Outliers  are  once  again 
observed  in  the  GDEM  data,  the  larger  out¬ 
lier  in  the  ASUW  scenario,  and  the  greater 
number  of  outliers  in  the  ASW  scenario.  The 
wide  dispersion  of  the  GDEM  derived  cov¬ 
erage  indicates  that  weapon  effectiveness  will 
vary  depending  on  location.  This  is  similar 
to  the  predictions  for  the  shallow  depth  band 
and  would  indicate  that  GDEM  predicts  a 
water  column  that  has  varying  coverage  val¬ 
ues  depending  on  horizontal  and  vertical 
location.  MODAS  data  once  again  indicates 
an  overall  performance  in  the  region  that  is 
slightly  less  than  the  GDEM  prediction; 
however,  the  MODAS  data  is  grouped  even 
more  tightly  than  in  the  shallow  depth  band. 
The  coverage  in  the  ASW  scenario  in  par¬ 
ticular  varies  little  about  the  mean  value.  This 
and  the  shallow  depth  band  predictions  in¬ 
dicate  uniform  coverage  can  be  expected 
even  at  some  depth. 

In  the  deep  depth  band  the  graphs  take 
on  a  slighdy  different  shape,  but  they  con¬ 


vey  much  the  same  meaning  (Fig.  12).  In 
both  scenarios  the  GDEM  graphs  are 
weighted  heavily  to  the  right  end,  predicting 
that  in  the  deep  depth  band  coverage  will  be 
very  good  over  most  of  the  area.  The  ASUW 


6.2.  Difference  of  IVIK48  Acoustic 
Presets  Using  GDEM  and  MODAS 

For  the  shallow  depth  band,  the  RMSD 
in  the  percentage  coverage  area  was  small 
over  most  of  the  AOI,  consistent  with  the 
similar  means  and  range  of  values  noted  in 
the  previous  section  (Fig.  13).  The  areas 


FIGURE  12 

Deep  depth  band  coverage  percentage  distributions:  upper  panels  for  GDEM  and  lower  panels  for  MODAS,  left 
panels  for  ASUW  scenario  and  right  panels  for  ASW  scenario. 
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FIGURE  14 

RMSD  for  medium  depth  band  coverage:  left  panel  for  ASUW  scenario  and  right  panel  for  ASW  scenario. 


computed  to  have  small  RMSD  coverage 
percentages  also  had  small  RMSD  in  tem¬ 
perature  and  salinity.  In  the  region  where 
the  RMSD  in  temperature  and  salinity  was 
largest,  though,  a  large  RMSD  in  percent¬ 
age  coverage  is  also  observed.  These  larger 
values  are  likely  areas  where  the  GDEM 
data  generates  overly  optimistic  coverage 
percentage  predictions.  For  the  surface  sce¬ 
nario,  RMSDs  of  up  to  25%  are  shown  in 
the  region  around  39°  N  73°  W,  and  the 
warm  salty  intrusion  observed  on  the 
MODAS  data  coincides  with  a  second  peak 
in  the  northeastern  section  of  the  graph. 
Overall  the  ASW  scenario  shows  RMSDs 
that  are  similar  to  the  ASUW  scenario,  the 
only  difference  being  that  the  values  are, 
on  average,  slightly  smaller.  The  notable 
exception  is  the  peak  located  at  the  top 
portion  of  the  graph. 

For  the  mid  depth  band  the  percentage 
coverage  RMS  difference  for  the  ASUW  sce¬ 
nario  is  simply  a  scaled  down  version  of  the 
shallow  depth  band  ASUW  graph  (Fig.  14). 
This  makes  a  great  deal  of  sense  considering 
the  fact  that  the  coverage  percentage  distri¬ 
butions  for  the  shallow  and  mid  depth 
ASUW  scenarios  were  very  similar.  The  real 
difference  is  in  the  ASW  scenario.  The  single 
exceptional  peak  at  the  top  of  the  previous 
graph  is  gone  and  the  observed  differences 
have  become  much  smaller.  Most  of  the 
RMS  differences  for  the  mid  depth  ASW 
scenario  do  not  exceed  10%.  This  is  prob¬ 
ably  due  to  nearly  identical  coverage  per¬ 


centage  means  from  both  data  sets,  and  the 
tighter  grouping  of  the  GDEM  data  cover¬ 
age  percentage  predictions  in  the  mid  depth 
band  ASW  scenario.  The  RMSD  values  are 
small  even  in  the  areas  where  the  tempera¬ 
ture  and  salinity  differences  were  observed 
to  be  large,  such  as  in  the  upper  section  of 
the  graph. 

The  RMSD  observed  in  the  deep  depth 
band  scenarios  are  smaller  than  those  of  the 
shallow  depth  band,  but  similar  in  magni¬ 
tude  to  the  mid  depth  band  (Fig.  15).  For 
the  ASUW  scenario  the  RMSD  peaks  near 
the  northwestern  corner  of  the  AOI  and  then 
decreases  steadily  in  steps  heading  toward 
the  opposite  corner.  While  the  individual 


RMSD  values  seen  are  not  as  large  as  some 
of  the  ones  in  the  other  depth  bands,  more 
of  the  area  has  a  non-negligible  RMSD.  The 
cause  of  this  can  be  seen  from  the  percent¬ 
age  coverage  distribution  for  the  deep  depth 
ASUW  scenario. 

The  GDEM  data  results  in  values  that 
are  almost  all  larger  than  the  largest 
MODAS  derived  values.  This  overly  opti¬ 
mistic  prediction  means  that  over  a  large 
portion  of  the  AOI,  the  RMSD  is  going  to 
be  non-zero.  The  RMS  difference  in  the 
ASW  scenario  changes  very  little  from  the 
mid  depth  band  save  for  the  fact  that  the 
values  in  the  lower  right  corner  are  smaller. 
The  coverage  distributions  for  the  deep 
ASW  scenario  are  similar  to  the  ASUW 
case,  but  the  separation  between  the  two 
means  is  not  so  pronounced.  The  result  is 
a  larger  region  where  the  RMSD  is  small 
or  zero. 

Both  of  these  graphs  match  the  pattern 
that  has  so  far  been  observed  in  the  other 
depth  bands.  The  ASUW  scenario  has  the 
higher  RMS  difference  values,  with  areas  of 
both  high  temperature  and  salinity  differ¬ 
ences  corresponding  to  peaks  on  the  graphs. 
The  RMS  difference  values  also  approach 
zero  moving  toward  the  top  left  or  bottom 
right  corners.  Also,  as  depth  increases,  the 
difference  between  the  two  data  sets  de¬ 
creases  causing  the  difference  between  the 
coverage  percentages  to  decrease. 


FIGURE  13 

RMSD  for  shallow  depth  band  coverage:  left  panel  for  ASUW  scenario  and  right  panel  for  ASW  scenario. 
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FIGURE  15 

RMSD  for  deep  depth  band  coverage:  left  panel  for  ASUW  scenario  and  right  panel  for  ASW  scenario. 


Conclusions 

By  looking  at  the  RMSD  in  the  tem¬ 
perature  and  salinity  fields  generated  from 
the  GDEM  and  MODAS  data,  it  is  pos¬ 
sible  to  look  for  areas  where  the  data  differ 
significantly.  It  is  at  these  points  that  the  dif¬ 
ference  in  the  preset  effectiveness  should  be 
the  greatest.  This  is  observed  for  both  sce¬ 
narios  at  all  depth  bands.  The  percentage 
coverage  is  the  most  different  at  points  where 
both  the  temperature  and  salinity  RMS  dif¬ 
ferences  are  large.  This  is  especially  true  for 
the  shallow  depth  band  where  differences 
of  25%  are  observed  for  both  scenarios.  It  is 
of  interest  to  note  that  even  at  the  surface 
the  RMS  differences  for  the  temperature  and 
salinity  are  never  more  than  a  few  degrees  or 
PSU.  Even  with  only  this  slight  increase  in 
the  accuracy  of  the  inputs,  a  large  increase 
in  the  accuracy  of  the  prediction  of  the 
weapon  effectiveness  occurred.  This  seems 
to  imply  that  the  sensitivity  of  the  presets  to 
changes  in  the  inputs  is  quite  high. 

From  the  output  distributions  it  becomes 
clear  that  the  GDEM  derived  coverage  per¬ 
centages  indicate  that  weapon  effectiveness 
should  vary  not  only  in  the  horizontal  but 
also  in  the  vertical.  The  implication  is  that 
in  some  areas  coverage  will  be  very  high  and 
in  others  the  coverage  will  be  very  poor,  but 
the  tendency  is  for  the  coverage  to  be  high 
for  any  given  area.  The  MODAS  derived 
percentages  reveal  that  the  exact  opposite  is 
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true.  The  coverage  will  be  consistent  no 
matter  what  the  horizontal  location  or  depth 
band.  This  is  an  important  result  since  pre¬ 
diction  of  weapon  effectiveness  is  vital  to 
mission  planning  and  execution.  In  this  case 
an  unrealistic  expectation  in  the  weapons 
effectiveness  would  have  resulted  from  the 
use  of  the  GDEM  data  to  predict  the  cover¬ 
age  percentages  in  the  water  column.  The 
MODAS  data  also  would  have  given  the  user 
the  freedom  to  operate  anywhere  in  the  re¬ 
gion  knowing  that  their  weapon  would  func¬ 
tion  about  the  same  no  matter  the  location. 

The  most  obvious  limitation  of  this  work 
was  the  limited  data  set.  Any  future  work 
should  include  data  that  covered  a  wider 
number  of  areas  and  times.  Areas  of  strong 
thermal  and  salinity  contrast  are  of  particu¬ 
lar  interest.  Various  combinations  of  the  user 
inputs  into  the  WAPP  should  also  be  stud¬ 
ied.  The  effects  of  variables  such  as  bottom 
type  and  position  (upslope/downslope)  need 
to  be  addressed.  Another  avenue  of  study  is 
the  determination  of  how  the  number  of 
altimeters  affects  the  accuracy  of  the  out¬ 
puts.  It  has  been  determined  that  the  pre¬ 
sets  are  sensitive  to  the  addition  of  satellite 
data.  However,  the  effect  of  the  number  of 
satellite  inputs  still  remains  to  be  determined. 
Once  this  is  done  an  optimal  number  of  al¬ 
timeters  can  be  determined  based  on  mini¬ 
mizing  cost  and  maximizing  preset  accuracy. 


Acknowledgments 

This  work  was  jointly  supported  by  the 
Space  and  Naval  Warfare  System  Command 
and  the  Naval  Postgraduate  School.  The 
authors  would  like  to  thank  Charlie  Barron 
at  the  Naval  Oceanographic  Office  for  pro¬ 
viding  MODAS  data  set. 


References 

Chu,  P.C.,  E.  Gottshall  and  T.E. 
Halwachs.1998.  Environmental  effects  on 
Naval  warfare  simulations.  Institute  of  Joint 
Warfare  Analysis.  Naval  Postgraduate  School. 
Technical  Report,  NPS-IJWA-98-006.  33  pp. 

Chu,  P.C.,  S.K.  Wells,  S.D.  Haeger,  C. 
Szczechowski  and  M.  Carron.  1 997.  A 
parametric  model  for  the  Yellow  Sea  thermal 
variability.  J  Geophys  Res.  102:5655-5668. 

Chu,  P.C.,  Q.Q.  Wang  and  R.H.  Bourke. 
1999.  A  geometric  model  for  Beaufort/ 
Chukchi  Sea  thermohaline  structure.  J  Atmos 
Oceanic  Technol.  16:613-632. 

Fox,  D.N.,  WJ  Teague,  C.N  Barron,  M.R. 
Carnes  and  C.M.  Lee.  2002.  The  Modular 
Ocean  Data  Assimilation  System  (MODAS). 

J  Atmos  Oceanic  Technol.  19:240-252. 

Gandin,  L.S.  1965.  Objective  analysis  of 
meteorological  fields.  Israel  Program  for 
Scientific  Translation.  242  pp. 

Gottshall,  E.  1997.  Environmental  Effects  on 
Naval  warfare  simulations.  MS  Thesis  in 
Physical  Oceanography,  Naval  Postgraduate 
School,  (Monterey).  92  pp. 

Jacobs,  G.A.,  C.N.  Barron,  M.R.  Carnes, 
D.N.  Fox,  H.E.  Hurlburt,  P.  Pistek,  R.C. 
Rhodes,  W.J.  Teague,  J.P  Blaha,  R.  Crout, 
O.M.  Smedstad  and  K.R.  Whitmer.  1999. 
Naval  Research  Laboratory  Report  NRL/FR/ 
7320-99-9696.  Navy  Altimeter  Data 
Requirements.  25  pp. 

Teague,  W.J.,  M.J.  Carron  and  P.J.  Hogan. 
1990.  A  Comparison  between  the  Generalized 
Digital  Environmental  Model  and  Levitus 
Climatologies.  J  Geophys  Res.  95:7167-7183. 


